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FOUNDATIONS FOR PROPER-TIME RELATIVISTIC 
QUANTUM THEORY 


T. L. GILL, T. MORRIS, AND S. K. KURTZ 

Abstract. This paper is a progress report on the foundations for the 
canonical proper-time approach to relativistic quantum theory. We first 
review the the standard square-root equation of relativistic quantum 
theory, followed by a review of the Dirac equation, providing new in¬ 
sights into the physical properties of both. We then introduce the canon¬ 
ical proper-time theory. For completeness, we give a brief outline of 
the canonical proper-time approach to electrodynamics and mechanics, 
and then introduce the canonical proper-time approach to relativistic 
quantum theory. This theory leads to three new relativistic wave equa¬ 
tions. In each case, the canonical generator of proper-time translations is 
strictly positive definite, so that each represents a true particle equation. 

We show that the canonical proper-time version of the Dirac equation 
for Hydrogen gives results that are consistently closer to the experimen¬ 
tal data, when compared to the Dirac equation. However, these results 
are not sufficient to account for either the Lamb shift or the anomalous 
magnetic moment. 
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J This paper is an extended version of the talk given at the The 9th Biennial Confer¬ 
ence on Classical and Quantum Relativistic Dynamics of Particles and Fields held at the 
University of Connecticut in June 2014. 
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Introduction 

Following Dirac’s quantization of the electromagnetic field in 1927, and 
his relativistic electron theory in 1928, the equations for quantum electrody¬ 
namics QED were developed by Heisenberg and Pauli in the years 1929-30. 
From the beginning, when researchers attempted to use the straightforward 
and physically intuitive time-dependent perturbation expansion to compute 
physical observables, a number of divergent expressions appeared. Although 
it was known that the same problems also existed in classical electrodynam¬ 
ics, Dirac had shown that, in this case, one could account for the problem of 
radiation reaction without directly dealing with the self-energy divergence by 
using both advanced and retarded fields and a particular limiting procedure. 
Early attempts to develop subtraction procedures for the divergent expres¬ 
sions were very discouraging because they depended on both the gauge and 
the Lorentz frame, making them appear ambiguous. These problems were 
solved via the fundamental work of Feynman, Schwinger, and Tomonaga. In 
recent times, it is generally agreed that quantum electrodynamics (QED) is 
an almost perfect theory, which is in excellent agreement with experiment. 

The fact that QED is very successful is without doubt. However, there are 
still some foundational and technical issues, which require clarification and 
which leaves the thoughtful student with a sense of unease in taking this as 
the final answer. In light of the tremendous historical success of eigenvalue 
analysis in physics and engineering, it is not inappropriate to reinvestigate 
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the foundations with an eye towards clearly identifying the physical and 
mathematical limitations to our understanding of the hydrogen spectrum as 
an eigenvalue problem. 

In the first section of this review we take a new look at the square-root 
operator and show that it has an analytic representation as a nonlocal com¬ 
posite of three singularities. The particle component has two negative parts 
and one (hard core) positive part, while the antiparticle component has two 
positive parts and one (hard core) negative part. This effect is confined 
within a Compton wavelength such that, at the point of singularity, they 
cancel each other providing a finite result. Furthermore, the operator looks 
like the identity outside a Compton wavelength. 

In the second section, we provide an analytic diagonalization of the Dirac 
operator. Our approach leads to a complete split of the particle and an¬ 
tiparticle parts into two non-hemitian components, which are mapped into 
each other by the charge conjugation transformation. Thus, the full matrix¬ 
valued operator is Hermitian and shows (as is explained in the text) that the 
spinor representation in the Dirac equation hides its time nonlocal property. 
We conclude that the Dirac and square-root operator do not represent the 
same physics, despite that fact that they are related by a unitary transfor¬ 
mation. 

In the third section, we introduce the canonical proper-time approach 
to electrodynamics and mechanics. By convention, this approach fixes the 
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proper-time of the observed system as the clock of choice for all observers 
and explicitly shows that the question of simultaneity is actually a question 
of clock conventions. The change in convention produces a new symme¬ 
try group which is distinct from, but closely related to the Lorentz group, 
but with a Euclidean representation space. Thus, the new convention also 
replaces the standard form of Lorentz covariance by a new one. The advan¬ 
tage is that, this allows us to construct a parallel image of the conventional 
Maxwell theory for a charged particle, which is mathematically, but not 
physically, equivalent to the conventional form. The new wave equation 
contains a gauge independent term, which appears instantaneously along 
the direction of motion, but opposing any applied force and is zero other¬ 
wise. This is the near field (i.e., the field at the site of the charged particle). 
This shows that the origin of radiation reaction is not the action of a charge 
on itself but arises from inertial resistance to changes in motion. We show 
that the dissipative term is equivalent to an effective mass so that classical 
radiation has both a massless and a massive part. We also discuss solutions 
to a number of other problems that are solved with our new convention, 
which are either impossible or problematic within the standard framework. 

In the forth section, we describe the canonical quantized proper-time the¬ 
ory. We obtain three possible relativistic wave equations, because of new 


possibilities, for the manner in which the potential energy may be introduced 
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into the theory. Each new equation is generated by a strictly positive defi¬ 
nite canonical Hamiltonian, so that it represents a consistent particle. We 
focus on the proper-time extension of the Dirac equation. A basic test of our 
proper-time theory is the extent that it compares to the Dirac theory in ac¬ 
counting for the hydrogen spectra. We show that our canonical proper-time 
version of the Dirac equation gives results which are consistently closer to 
the experimental data, when compared to the Dirac equation. The present 
theory has not yet accounted for the Lamb shift or the anomalous magnetic 
moment. However, the analysis in sections one, two and three support our 
contention that the electron is not a point particle. This non-point nature is 
only expected to be important in s-states, where there is a finite probability 
of the electron being at the center of the proton. This aspect of our research 
is still in progress and will be reported on at a later time. 


1. The square-root equation 

In the transition to relativistic mechanics, the equation E 2 = c 2 p 2 + m 2 c 4 
leads to the quantum Hamiltonian 

H = y/ c 2 p 2 + m 2 c 4 . 


Thus, it is quite natural to expect that the first choice for a relativistic wave 
equation would be: 


,h m = 


\J C 2 [p — i./VlA 1 I tu-r‘ + V 


V 7 , 
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where p = —ihX7. However, no one knew how to directly relate this equation 
to physically important problems. Furthermore, this equation is nonlocal, 
meaning, in the terminology of the times (1920-30), that it is represented 
by a power series in the momentum operator. One was led in this way to 
the Gordon-Klein and Dirac equations. 

1.1. Background. Since the early work, many investigators have studied 
the square-root equation. It is not our intention to provide a detailed his¬ 
tory or to identify the many important contributors to the study of this 
problem. In recent times, the works of Silenko (see m, m are well worth 
reading. They also provide a very good list of the important historical stud¬ 
ies. In addition, he has made a number of interesting investigations into the 
transformational relationship between the square-root and Dirac equation 
m is a good starting point). The recent paper by Simulik and Krivsky 
g] offers another interesting approach to the square-root equation and its 
relationship to that of Dirac. Closer to our investigation of the square-root 
equation is the study by Kowalski and Rembielinski [5] (also known as the 
Salpeter equation). They have used it as an alternative of the Klein-Gordon 
equation. 

In this section, we take a new look at the square-root equation. First, 
we investigate the extent that the non-commutativity of p and A affect 
our ability to give an unambiguous meaning to the square-root operator. 
We show that a unique analytic representation is well defined for suitable 
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time-independent A provided we can solve a corresponding equation of the 
Schrodinger type. We then investigate a few simple cases of solvable models 
in order to get a feeling for the physical interpretation of this operator. 

To begin, we start with the equation: 

(1.1) S 1 ^] = Hs'tP = |/3y / c 2 (p — f A) 2 - ehcT, • B + m 2 c 4 | ip. 

Where /3 and £ are the Dirac matrices 



I 0 


1- 

O 

b 

i_ 

II 

1 

O 

1 

hH 

, £ = 

0 a 


I and cr are the identity and Pauli matrices respectively. Under physically 
reasonable mathematical conditions, the following operator is a well defined 
self-adjoint generator of a strongly continuous unitary group: 

H 2 = c 2 [p — (e/c)A] 2 — ehcTj ■ B + m 2 c 4 . 

From the basic theory of fractional powers of closed linear operators, it can 
be shown that 

(1.2) v / ^| = (v/fl? y 1 H? = Hl(yH?)~ 1 . 

In order to construct an analytic representation for equation (1.1), we as¬ 
sume that B is constant. The general case can be found in [6]. Let 
G = — c 2 (p — f A) 2 and uj 2 = m 2 c 4 — ehc£ ■ B, so that uj is also constant. 
Using this notation, we can write (1.1) as 


= |W-G + w 2 } ^ 
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Using the analytic theory of fractional powers of closed linear operators and 
equation (1.2), we can represent Sty] as 

a r °° i d\ 

(1.3) Sty] = — [(A + u 2 ) — G] (—G + u 2 )—=ty], 

71 Jo V A 

where [(A + w 2 ) — G] 1 is the resolvent associated with the operator 
—G + cj 2 . The resolvent can be computed directly if we can find the funda¬ 
mental solution to the equation: 

d Q(x, y; t)/dt + (G - uj 2 )Q(x, y; t) = 5(x - y). 


Schulman [7J has shown that the solution to the above equation is 


/• x (t) f rt r x 

Q= ^[x(s)] exp / U[x(s)]ds + jg / A[x(s)] • dx(s) 
Jx(0) l ■/ 0 Jy 


where V = /ifr and 


/•x(t)= x />x(t)=x f rt 

/ Z>[x(s)] = / ^[ x (s)] exp < / 

lx(0)=y J x(0)=y [ J 0 


dx(s) 


ds 


ds 


= lim 

N—too 


■j nJV /2 


N 


N 


47 T 


II dxj exp <j - [ 4 ^ (Xj - Xj^f 

k =1 | j =1 


and £n = t/N. A rigorous justification for the path integral can be found 
in Gill and Zachary [ 8 ]. We assume that f* A[x(s)] • dx(s) = A • (x - y), 
where A is the mean value of A. Using this, we have: 


/■= 

[(A + J 2 ) - G] f(x) = 

Jo 


—A t 


Q(x,i;y,0)/(y)dy 


dt 
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and 


[(A H- co 2 ) — G] /(x) 

= [ e {£ A <*-U 

J K 


exp 


(x-y)" _ io 2 t _ \ f 

At ~W AL 


dt 


J (Airt) 


Tj2 \ f{y)dy 


Using a table of Laplace transforms, the inner integral can be computed to 
get 

i - ~ -i 

dt 


oo 

exp 

(x - y ) 2 w 2 t 

4 t h 2 Xt 

1 exp 

-v / ( a + m 2 ) x - y 


(4vrt) 


3/2 


4vr ||x — y || 

where pi 2 = u 2 /h 2 . Equation (1.3) now becomes 


^M(x) 


Once again, we interchange the order of integration and perform the com¬ 
putations to get 


poo [ exp 

-V / ( a + M 2 ) l x - y || 

1 d\ 4/rr(|)K 1 [pi x — y ] 

/o { 

x - y 

J VA 7T 1 / 2 ||x — y || 


where Ki [z] is the modified Bessel function of the third kind and first order. 
Thus, if we set, a = and a = we get 

SM(x) = $ f e^-y)] ^H x ~ y IU (—G + u 2 )^y)dy 

(1.4) ^ l|X ' yl1 

= fy~G +UJ 2 ) [ e^'^-y)} ^ Kll|X ~ y 1 II] ^(y)rfy. 

Jr 3 ll x -y|| 

Since V • a = 0, we have 


—G + J 2 = h 2 (—A + 2ia • V + a 2 + pi 2 ) , 
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so that so that (1.4) becomes 

SM(X) = 4^ (-A + 2 za • V + a 2 + p 2 ) J * e -(x-y)MKifcdO ^ (y)dy . 

The operator (—A + 2ia • V + a 2 + /x 2 ) acts on x, making the integral sin¬ 
gular. However, this singular representation constructed below, equation 
(1.5) has many of the properties observed in experiments. As will be seen, 
it represents the confinement of three singularities within a Compton wave¬ 
length. (A full discussion is delayed to the end of this section.) 

We omit many of the computational details, which can be found in [9], 
but the idea is to consider a ball B p (x) of radius p about x, so that M 3 = 
M 3 U Bp(x), where M 3 = (M 3 \B p (x)) and tM 3 = (cM 3 \<9Bp(x)). We then 
restrict all operations to M 3 and only let p —>• 0 at the end. 


1.2. Free Case. The free particle case is the simplest (but still interesting), 
with A = 0, so that 


5[^](x) 

(1.5) 

_ ffih 2 c/3 

— 


/ 

R 3 


l|x-y|| 


— 47 rh (x — 


< x-y )] 


Ml x ~y||] I 2Ki 0||x—y|| 


l|x-y|| 


+ 


Hl x -y|| 


■} ^(y)rfy- 


If x 7 ^ y, the effective kernel of equation (1.5) is 


Ao[/x||x-y||] , 2Ki [p ||x — y| 

,, O “T 


I x -yir Hl x -y| 
1-2 


Recall that the integral of ||x — y|| over M 3 is finite. In order to understand 
the physical interpretation of equation (1.5), it will be helpful to review 


some properties of the modified Bessel functions Ko[u], u 2 Ad/ 2 M and 
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u 1 K\[u]. We follow Gradshteyn and Ryzhik jl()j . for 0 < u <C 1, we have 
that: 

= [l + 0i (u)]u~ 2 

A l/2 M _ - 1 

u l/2 ~\J 2 U 

Kq [u] = [1 + 0oO)] lnw _1 , 

where Qq, 0\ —>• 0 as n —>• 0. We note that, up to a multiplicative constant, 
_i 

u 2 K ] / 2 [u] is the well-known Yukawa potential [TTJ, conjectured in 1935 
to account for the short range of the nuclear interaction. From here, we 
see that, near u = 0, the singular term u~ 1 K\[u] is twice as strong as the 
Yukawa potential. The singular term Kq[u] is actually integrable and so 
does not contribute at u = 0. Looking at equation (1.5), we see that the 
singular term —87 xKi{u)5(u) acts to cancel the the singular term u -1 /\i[u] 
at u = 0, so that the total integral is well defined. 

The behavior of these functions is quite different, when u> 1. In this 
case, we have: 

[ 1 + ^ 

-fc'i/2 M _ exp {-it} 

u l/2 V 2 y, 

A-„ [„] = [1 + «i(«)] 

We see that each term has a exponential cutoff. However, now Kq[u] has the 
longest range, while u^K^u] has the shortest range. Furthermore, inspec¬ 
tion shows that is multiplied by the reduced Compton wavelength, 


Ki M 

u 


which further shortens its range. 
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It is clear that our square-root operator represents an extended object 
with an effective extent of about a Compton wavelength. 


1.3. Constant A Case. When A / 0 is constant, V • A = 0 and B = 
V x A = 0, so that we get (a = a): 


S[V>](x) = 


-p?T?cP 

-^ 


l f 

R 3 


(x-y) 


ll x —y II 


>-y)] { 


1 _ 4^(5 ( x _ y ) \ K o[H|x—yll ] _|_ 2KiM|x-y|| ] 


l|x—y II 


H|x-yr 


} V>(y )dy- 


1.4. The Constant Field Case. If B ^ 0 is constant, then A(z) = ^zxB. 
Let a(z) = ^|-z x B and F = —a • (x — y). In this case A(z) • dz = 0, so we 
can write the final result as: 




/R 3 


1 _ 47r5(x—y) 

||x—y || 1+iF 


[1 + iF] Sgtf 


( 16 ) 

+ £gl a 2 Kl g7» lV(y)rfy, 

J R 3 


where 


^ ^ K 2 [/x ||x — y||] = K 0 [/i ||x — y ||] + 2K X Hx-yjj] 

ll x _ y II ll x -y|l mII x - y || 2 

From equations (1.6) and (1.7), we see that a constant magnetic field makes 
a real difference compared to either the A = 0 or A / 0 cases, producing 
two extra terms, in addition to the free particle term. The first new term is 
purely imaginary and singular at x = y (like the Yukawa term). Physically, 
we interpret this term as representing particle absorption and emission (see 
Mott and Massey m . The second term is real, repulsive and nonsingular. 
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In addition, the effective mass fi is constant but matrix-valued with complex 
components, /i 2 = m 2 c 2 /h 2 — ^X • B. Since 



( 

\ 


f 

\ 


f 

\ 


( 

\ 


a 

0 


0 

l 


0 

—i 


1 

0 

X = 



; o-i = 



, 02 = 



, 03 = 




1° 

a ) 


l 1 

°J 


V * 

0 J 


l 0 

-v 


(T-/l S 3)l2 !(£ 2 -^i)I 2 

=^(B 2 -iB 1 )I 2 (^ + £B 3 ) I 2 


From known properties of Bessel functions for nonintegral u, we can repre¬ 
sent K„[it] as 

, 0/ I-„(«) - I„(«) _ e*/ 2 (^) J_„(*u) - e-'/ 2 ^) J„{iu) 

(2/7t).K„|mJ — . — 

sm 7r/y sm tiv 


In the limit as v approaches an integer, the above takes the indeterminate 
form 0/0 and is defined via L’Hopital’s rule. However, for our purposes, 
we assume that v is close to an integer and u = u\ + iu 2 , u 2 / 0. In this 
case, K v [it] acquires some of the oscillatory behavior of 3 v [u}. Thus, we can 
interpret equation (1.6) as representing a pulsating mass (extended object 
of variable mass) with mean value 7i/c||/i||, where ||/t|| = [/T/x] 1 / 2 and n* is 
the Hermitian conjugate of //, with the square root being computed using 
elementary spectral theory. If B is very large, we see that the effective mass 
can also be large. However, the operator still looks (almost) like the identity 
outside a Compton wavelength. 

In closing, we should say a few additional words about the interesting 
work of Kowalski and Rembielinski [5j. They solve the Salpeter equation 
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(/? = I) and construct a number of examples. This work represents an 
original contribution to our understanding of the square-root equation. They 
approach the problem using the method of Fourier transforms and get the 
correct solution for x / y, However, this approach misses the the x = y 
term, giving the impression that their equation is not defined at that point. 
This minor defect can be easily fixed, by adding our delta term, which 
makes their solutions well-defined for all x. More important, is to note 
that replacing the indentity operator by the (general) (3 matrix provides a 
generalization of their solutions for all spin-values. 


1.5. Conclusions. From our analysis, we have the following conclusions 
concerning the square-root operator: 

(1) In the simplest case, A = 0, the square-root operator has a repre¬ 
sentation as a nonlocal composite of three singularities. The particle 
component has two negative parts and one (hard core) positive part, 
while the antiparticle component has two positive parts and one 
(hard core) negative part. This effect is confined within a Compton 
wavelength such that, at the point of singularity, they cancel each 
other providing a finite result. Furthermore, the operator looks like 
the identity outside a Compton wavelength. (Recall that the ex¬ 
perimental observation of three singularities in proton and neutron 
scattering experiments led to the quark model.) 
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(2) A constant magnetic field induces changes in both the mass and 
the shape of this extended object. It also increases the number of 
singularities. This suggests that the square-root operator represents 
a charge/mass density, for otherwise it could not be affected by a 
constant magnetic field. 

(3) The square-root operator is not physically the same as the Dirac 
operator despite the fact that they are related by a unitary trans¬ 
formation. (We will discuss this point further in the next section.) 

2. The Dirac Equation 

The first successful attempt to resolve the question of how best to handle 
the square-root equation was made by Dirac in 1926 m Dirac noticed that 
the Pauli matrices could be used to write c 2 p 2 + m 2 c 4 as [ca • p + mc 2 /3] 2 . 
The matrix a is defined by a = (ai, 0 : 2 , 0 : 3 ), where 



/ 

\ 


( 



( 

\ 


f 

\ 


0 



0 

1 


0 

—* 


1 

0 

OLi = 



, cn = 



, CQ = 



, ^3 = 





0 ) 


l 1 

°J 


V i 

0 j 


1 ° 

- 1 / 


Thus, Dirac m showed that an alternative representation of the square-root 
equation could be taken as: 

(2.1) ^~Qt = ZOi ’ P mc2 ft\ 

In this case, T must be viewed as a vector-valued function or spinor. To be 
more precise, f 6 L 2 (R 3 ,C 4 ) = L 2 (R 3 ) <S> C 4 is a four-component column 
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vector T = (V’l, 4>2, </h, ■ In this representation, ?/> = (V’l,^2)* repre¬ 

sents the particle (positive energy) component, and tp = (<^ 1 ,^ 2 ) represents 
the antiparticle (negative energy) component of the theory (for details, see 
Thaller [14]). 

A fair understanding of the Dirac equation can only be claimed in recent 
times, and, as pointed out by D. Finkelstein, “Dirac introduced a Lorentz- 
invariant Clifford algebra into the complex algebra of observables of the 
electron”. (See, in particular, Biedenharn [15] or deVries [16] and Hestenes 
Em) Despite successes, both practical and theoretical, there still remain 
a number of conceptual, interpretational, and technical misunderstandings 
about this equation. It is generally believed that it is not possible to separate 
the particle and antiparticle components directly without approximations 
(when interactions are present). The various approximations found in the 
literature may have led to this belief. In addition, the historically important 
algebraic approaches of Foldy-Wouthuysen [18], Pauli m, and Feynman 
and Gell-Mann [20] have no doubt further supported such ideas. 

In this section we show that it is possible to directly separate the particle 
and antiparticle components of the Dirac equation without approximations, 
even when scalar and vector potentials of quite general character are present 
(see 0)- We show that the square root operator cannot be considered 
physically equivalent to the Dirac operator from another point of view. In 
addition, we offer another interpretation of the zitterbewegung and the fact 
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that the expected value of a velocity measurement of a Dirac particle at any 
instant of time is ±c. 


2.1. Complete Separation. It turns out that a direct analytic separation 
is actually quite simple and provides additional insight into the particle 
and antiparticle components. In order to see this, let A(x, t) and C(x) be 
given vector and scalar potentials and, after adding V (x) and making the 
transformation p — )• it = p — e/cA, write (2.1) in two-component form as: 
d'lb 

ih~— = (V + mc 2 )i/) + c(cr • 7 j)ip 

(2 ' 2) * 

ih~— = (V — mc 2 )(p + c(a ■ 
at, 

We write the second equation as: 

ip = —iDipi B\ = j t {V — m,c 2 ) and D = \c(cr ■ n). 

In this form, we see that from an analytical point of view equation (2.2) is 
a first order inhomogeneous partial differential equation. This equation can 
be solved via the Green’s function method if we first solve 


d 

m +iBl 


d_ 

m 


+ iB\ 


u(t ) 


<*(*)■ 


It is easy to see that the solution to this equation is 


so that 


u{t) = 6(t) exp{— iBit}, 9(t) 


1, t > 0 
0, t < 0 



cexp{— iB\{t — t)} [(ct • -k)/ ih] 
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It now follows via convolution that: 


(2.3) 


ih^j- = (V -\- me 2 )il> 
at 

+ [c 2 (cr ■ it) / ih\ f exp{—i£?i(t — r)}(<7 • 7r)^(r)dT. 


In a similar manner, we obtain the complete equation for cp: 

= (V - mc 2 )ip 

( 2 - 4 ) f t 

+ [c 2 (cr • 7T )/ih] / exp{— iB 2 (t — t)}(cj • 7 t)<p(t)cIt, 


where 


v(t) = 9(t.) exp{— iB 2 t}, 0(t) 


1, t > 0 
0, t < 0 


B 2 = \{v + me 2 ). 

Thus, we have decomposed L 2 (M 3 ,C 4 ) as L 2 (M 3 ,C 4 ) = L 2 (M 3 ,C 2 ) © 
L 2 (M 3 ,C 2 ). One copy of L 2 (M 3 ,C 2 ) contains the particle (positive energy) 
wave component, while the other copy contains the antiparticle (negative 
energy) wave component. Which of these copies corresponds to the compo¬ 
nents ^ = (VhjV^/ and which to the components p = (p 1 ,^ 2 )* depends, 
to some extent, on the properties of the scalar potential V. It may have 
been noticed that equations (2.3) and (2.4) are non-hermitian. It is shown 
in [9] that they are mapped into each other by the charge conjugation op¬ 
erator, so that the full matrix representation is hermitian. An unsettled 
issue is the definition of the appropriate inner product for the two sub¬ 
spaces, which will account for the quantum constraint that the total prob¬ 


ability integral is normalized. We can satisfy this requirement if we set 
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O, x) = ipiXi + ^ 2 X 2 , (ip, x)i = {Al'ip, Ax) and (ip, rj) 2 = (A 2 <p, A 2 rj), 
where Aiip = cu(t) * [(cr • ir)/ih] ip(t) and A 2 p> = cv(t) * [(a ■ ir)/ih] c p(t ). We 
can now define the particle and antiparticle inner products by: 


(2.5) 


(V> ixl= / [W>, x) + (V’, x)i] dx 

JR3 

<¥> i » 7 )ap = [(¥>> ? ?) + (¥>> ^7)2! dx > 


so that the normalized probability densities satisfy: 

rt 


( 2 . 6 ) 


pip = IV’I + 


/V = l^l 2 + 


/ cexp{— iB(t — t)} [(cr • 7r )/ih] ij)(T)di 

-OO 


’ —OO 

rt 


/ cexp{—i,B' (t — r)} [(cr • n)/ih] ip(r)dT 

-OO 


2.1.1. Interpretations. Writing the Dirac equation and the direct separation 
in two-component matrix form, we have: 


d 



(V + me 2 ) c(a ■ 7r) 



ih— 


— 




dt 







<p 


c(cr - it) (V — me 2 ) 


<p 


and 



(V + me 2 ) 


r 1 


0 


j- 



+ [c 2 (a ■ tt)/ ih] [u* (a ■ it)] 





(1 V — me 2 ) 


V 

L J 


0 




+ [c 2 (ct • n) j ih] [v * (a ■ 7r)] 



We call the latter equation the analytic diagonalization of the Dirac equation 
because the wave function has not changed. 


The standard approach to the diagonalization of the Dirac equation (with¬ 
out an external potential V) is via the Foldy-Wouthuysen representation 
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[18] . Assuming that A does not depend on t, the following generalization 
can be found in deVries M- 


■h 9 

lh m 


where 


$1 

d>2 


0 — ^/^Tr^-^cK{T^~By+rn?? 


S = 


l. ^ 


\0 <T ) 


In this case, 


t 


$1 $ 2 


Ufw 


*!> 


t 

and our square-root operator 


5 = u FW u n u^ v . 

From equation (2.2), we conclude that the coupling of the particle and 
antiparticle wave functions in the first-order form of the Dirac equation hides 
the second order nonlocal time nature of the equation. We know that the 
square-root operator is nonlocal in space. Thus, the implicit time nonlocality 
of the Dirac equation is mapped into the explicit spatial nonlocality of the 
square-root equation by the Foldy-Wouthuysen transformation. This is a 
mathematical relationship, which is not physically equivalent. 

The time nonlocal behavior raises questions about the zitterbewegung. 
The physically reasonable interpretation of the zitterbewegung and the fact 
that the expected value of a velocity measurement (of a Dirac particle) at 
any instant in time ±c are reflections of the fact that the Dirac equation 
makes a spatially extended particle appear as a point in the present by 
forcing it to oscillate between the past and future at speed ±c. 
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3. Classical Proper-time Theory 

In this section, we briefly review the classical theory. The theory was first 
introduced in 2001 [21] and further discussed in [22] . However, the theory 
has its roots in the foundations of quantum electrodynamics as developed 
by Feynman and Dyson. 


3.1. Background. Following the suggestions of Feynman and Dyson, our 
program began with the development of a mathematical theory for Feyn¬ 
man’s time-ordered operator calculus, where time is accorded its natural 
role as the director of physical processes. Briefly, our theory is construc¬ 
tive in that operators acting at different times actually commute (in the 
mathematical sense). This approach allows us to develop a general pertur¬ 
bation theory for all theories generated by unitary evolutions. We are also 
able to reformulate our theory as a physically motivated sum over paths as 
suggested by Feynman. Our purpose was to prove the last two remaining 
conjectures of Dyson concerning the mathematical foundations for QED (see 
[23]). (A. Salam confirmed Dyson’s first conjecture [24] . while S. Weinberg 
[25] confirmed his second one.) In particular, we showed that: 

(1) The renormalized perturbation series of quantum electrodynamics is 
at most asymptotic. (We also provided the remainder so that, in the 
mathematical sense, the expansion is exact.) 
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(2) The ultraviolet divergence of quantum electrodynamics is caused by 
a violation of the time-energy uncertainty relations (at each point in 
time). 

As a special case, our approach also provided the first rigorous mathematical 
foundation for the Feynman path integral formulation of quantum mechanics 
(see [8]). 

In the Feynman world-view the universe is a three-dimensional motion 
picture in which more and more of the future appears as time evolves. Time 
is a physically defined variable with properties distinct from those of the 
three spacial variables. This view is inconsistent with the Minkowski world¬ 
view, in which time is an additional coordinate for space-time geometry. 

With this inconsistency in mind, we began to investigate the possibility 
that an alternative formulation of both classical and quantum theory could 
exist, which encodes the Feynman world-view. We discovered the canonical 
proper-time approach to classical electrodynamics, in which the proper-time 
of the observed system is used as opposed to the proper-time of the observer. 


3.2. Maxwell’s equations. For the local-time version of Maxwell’s equa¬ 
tions, it is convenient to start with the standard definition of proper-time: 


(3.1) 


dr 2 = dt 2 — \dx 2 = dt 2 
c z 


1 - 


W" 


w = 


dx 

dt 


Motivated by geometry and the philosophy of the times, Minkowski sug¬ 


gested that we use the proper-time to define a metric for the space-time 
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implementation of the special theory of relativity. Physically, it is well- 
known that dr is not an exact one-form because a particle can traverse 
many different paths (in space) during any given r interval. This reflects 
the fact that the distance traveled in a given t interval depends on the forces 
acting on the particle. This also implies that the clock of the source carries 
additional physical information about the acting forces. In order to see this, 
rewrite equation (3.1) as: 

. , , <> l,o , o [ u 2 l <ix 

(3.2) dt 2 = dr 2 +—?dx~ = dr 2 1-| —- , u =— . 

cr [_ c J «t 

For any other observer, we have: 

(3.3) dt' 2 = dr 2 + —dx 2 = dr 2 1 H- tt , u' = 

& |_ c \ “ r 

It follows that observers can use one unique clock to discuss all events as¬ 
sociated with the source (simultaneity). We also note that, the phase space 
variables remain unchanged because the momentum p = mw = mou, where 
rn = ymo- 

From equations (3.2) and (3.3) we see explicitly that, the new metric for 
each observer is exact, while the representation space is now Euclidean. In 
order to clearly see that we have a change in the clock convention, assume 
that we are observing a particle moving with constant velocity relative to 
our unprimed (inertial) frame. In this case, we can integrate equation (3.2) 
obtaining: 
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The inverse relationship is 



where w = Now, t and r differ by a scale factor, so that either may be 
used (a convention). The advantage of the r representation is that the same 
r is also available to our prime observer in her frame: 



The important point of our theory is that, this convention is available in 
the non constant velocity case. (In order to show that general case does not 
complicate matters, in Section 3.3 we construct the transformation group 
for all cases.) 

In our new formalism, the natural definition of velocity is no longer w = 
dx/df but u = dx/dr. This suggests that there may be a certain duality 
in the relationship between t. r and w, u. To see that this is indeed 
the case, recall that u = w J 1 — (w 2 / c 2 ). Solving for w, we get that 
w = u j y^l + (u 2 /c 2 ). If we set b = + u 2 , this relationship can be 

written as 


w 

c 


u 

b‘ 


(3.4) 
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For reasons to be clear momentarily, we call b the collaborative speed of 
light. Indeed, we see that 

Id _ 1 dr d _ 1 1 d _ 1 d 

cdt c dt dr c _|_ ( u 2 j c 2 ) 9t b dr' 

For our prime observer, it is easy to see that the corresponding result will 
be: 


nn — = — IJL = IJL 

1 J c ir Cdt' V dr 

From equations (3.5) (and (3.6)) we see that the non-invariance of t, {t') and 
the invariance of c on the left is replaced by the non-invariance of b,(b') and 
the invariance of r on the right. These equations represent mathematically 
equivalent relations. Thus, wherever they are used consistently as replace¬ 
ments for each other, they can’t change the mathematical relationships. In 
order to see their impact on Maxwell’s equations, in c.g.s. units, we have: 


V • B = 0, 


V • E = 47rp, 


(3.7) 


„ „ 1 dB „ „ 1 

V x E =- — , V x B = - 

c dt c 


3E 

— + 4»pw 


Using equations (3.1) and (3.2) in (3.3), we have the mathematically identical 
representation for Maxwell’s equations: 


V • B = 0, V • E = 47 rp, 

dE , 

—-1- 47T/5U . 

dr 

Thus, Maxwell’s equations are equally valid when the local time of the par¬ 
ticle is used to describe the fields. This leads to the following conclusions: 


V x E = - 


IdB 

blh 1 


V x B = 
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(1) There are two distinct clocks to use in the representation of 
Maxwell’s equations. (The choice of clocks is a convention.) 

(2) Since the two representations are mathematically equivalent, math¬ 
ematical equivalence is not the same as physical equivalence. 

(3) When the proper-time is used, the constant speed of light c is re¬ 
placed by the effective speed of light b, which depends on the motion 
of the system (i.e., b = \J c 2 + u 2 ). Thus, we have a natural varying 
speed of light theory (VSL), as opposed to a postulated one ( see 
Magueijo [26] or Moffat [27]). 


Let us now derive the corresponding wave equations in the local-time vari¬ 
able. Taking the curl of the last two equations Maxwell equations (above), 
and using standard vector identities, we get: 


1 5 2 B u • a 


(3.8) 


b 2 dr 2 


6 4 


19 2 E u • 


b 2 5r 2 


b 4 


5B 

77 

5E' 

57 


- V 2 • B = - [4vrV x (pu)\ , 
b 

-V 2 -E = -V(4np)-~ 


47t(/9u) 


where a = du/dr is the effective acceleration. The new (gauge independent) 
term appears instantaneously along the direction of motion, but opposing 
any applied force and is zero otherwise. This is the near field (i.e., the field 
at the site of the charged particle). This is exactly what one expects of the 
back reaction caused by the inertial resistance of a particle to accelerated 
motion and, according to Wheeler and Feynman [28], is precisely what is 
meant by radiation reaction. Thus, the collaborative use of the observer’s co¬ 


ordinate system and the local clock of the observed system provides intrinsic 
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information about the local field dynamics not available in the conventional 
formulation of Maxwell’s theory. It is shown in [22j, that the theory does 
not require point particles, self-energy divergence, mass renormalization or 
the Lorentz Dirac equation. 

It is also shown in [2Q that, for a closed system of interacting charged 
particles, the proper-time of the center of mass corresponds to the historical 
clock of Horwitz, Piron, and Fanchi (see [29] and [30]). In this case, b = c and 
the corresponding Maxwell equations represent the far field (only retarded 
potentials). It was further shown that, from this vantage point, the particle 
interactions appear as the delayed action-at-a-distance type. This verifies 
the Wheeler-Feynman conjecture that field theory and delayed action-at-a- 
distance are complimentary manifestations of the same physics (see [28]). 
The requirement of total conservation of momentum, angular momentum 
and energy allowed us to prove the complete absorption of radiation by all 
particles in the system. (Recall that, this was an assumption in the Wheeler 
and Feynman approach and the center-fold of their theory .) 

If we make a scale transformation (at fixed position) with E —>• (6/c) 1 / 2 E 
and B —>• (b/c) 1 ^ 2 B, the equations in (3.4) transform to 


1 a 2 B 

B + 

b 

3 b 2 

b 2 dr 2 

2 b 3 

4 b A 

(3.9) 


L 


1 <9 2 E ^ 

E + 

b 

3 b 2 

b 2 dr 2 

2 b 3 

4 6 4 


C V 2 

B= tf /2 t 47rV x ()1 ’ 


E = - 


M/2 

—j- 7 ^ V (47Tp) 


C^d_ 

6 3 / 2 dr 


47t(/9u) 

b 
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This is the Klein-Gordon equation with an effective mass // given by 
h 2 b 3 b 2 | f h 2 u • ii + u 2 5 (u • u) 2 | 

^2 J = Y# 2 ¥ IP j ' 

Remark 3.1. IFe note that, when b is constant, a = 0, fi = 0 and t = } jj 
(also t! = --t), so that the local time theory is both mathematically and 
physically equivalent to the standard theory. However, when b is not constant 
fi / 0 and the two approaches are not physically equivalent. 


For additional insight, let (x(t),t) represent the held position and 
(x(r / ),r / ) the retarded position of a source charge e, with r = x — x. If 
we set r = |x — x|, s = r — ( ^ r fe u ' ) ), and r u = r — fu, then we were able to 
compute the E and B fields directly in [21] to obtain: 

e[r u (l-u 2 /b 2 )] e [r x (r u x a)] e(u • a) [r x (u x r)] 

E(x - T) = -?- + —P?— + -- 

and 

e [(r x r u )(l - u 2 /& 2 )] er x [r x (r u x a)] er(u • a)(r x u) 
,T rs 3 r6 2 s 3 fe 4 s 3 


(It is easy to see that B is orthogonal to E.) The first two terms in the 
above equations are standard, in the (x(t),w(t)) variables. The third part 
of both equations is new and arises because of the dissipative term in our 
wave equation. (Once again, this term is zero when b is constant.) It is 
easy to see that r x (u x r) = r 2 u — (u • r)r, so we get a component along 
the direction of motion. (Thus, the E held has a longitudinal part.) This 
conhrms our claim that the new dissipative term is equivalent to an effective 
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mass that arises due to the collaborative acceleration of the particle. This 
means that the cause for radiation reaction comes directly from the use of 
the local clock to formulate Maxwell’s equations. It follows that, in this 
approach there is no need to assume advanced potentials, self-interaction, 
mass renormalization and the Lorentz-Dirac equation in order to account 
for it (radiation reaction), as is required when the observer clock is used. 
Furthermore, no assumptions about the structure of the source are needed. 

3.3. Proper-time Lorentz Group. We now identify the new transfor¬ 
mation group that preserves the first postulate of the special theory. The 
standard (Lorentz) time transformations between two inertial observers can 
be written as 

(3.11) t' = 7 (v) [t - x • v/c 2 ] , t = j(v) [t' + x' • v/c 2 ] . 

We want to replace t, t' by r. To do this, use the relationship between dt 
and dr to get: 

(3.12) t = l ( b(s)ds = ifer, t' = \ [ b'(s)ds = \b'r , 

Jo Jo 

where we have used the mean value theorem of calculus to obtain the end re¬ 
sult, so that both b and b' represent an earlier r-value of b and b' respectively. 
Note that, as b and b' depend on r, the transformations (3.8) represent ex¬ 
plicit nonlinear relationships between t,t r and r (during interaction). This 
is to be expected in the general case when the system is acted on by external 
forces. However, if b is constant (so is b'), then t, t' and r differ by a scale 
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transformation, which means they are physically equivalent, in addition to 
their natural mathematical equivalence. 

If we set 

d* = d/y(v) - (1 - 7 (v)) [(v • d)/( 7 (v)v 2 )] v, 

we can write the transformations that fix r as: 

x' = 7 (v) [x* — (v/c)6r] , x = y(v) [x 7 * + (v/c)6'r] , 

(3.13) u' = 7 (v) [u* - (y/c)b\ , u = 7 (v) [u* + (v/c)6'] , 

a' = 7 (v) {a* - v [(u • a)/(6c)]} , a = 7 (v) {a'* + v [(u 7 • a')/ (b'c)] } . 

If we put equation (3.8) in (3.7), differentiate with respect to r and cancel 
the extra factor of c, we get the transformations between b and b'\ 

(3.14) b'( t) = j(v) [6(r) — u • v/c], b(r) = j(v) [& 7 (r) + u 7 • v/c] . 

From these results, it follows that, at the local level, during interaction, 
equations (3.9) and (3.10) provide a nonlinear and nonlocal representation 
of the Lorentz group. We call it the proper-time Lorentz group. 

3.4. Proper-Time Particle Theory. We now investigate the correspond¬ 
ing particle theory. The key concept to our approach may be seen by exam¬ 
ining the time evolution of a dynamical parameter W (x, p), via the standard 
formulation of classical mechanics, described in terms of the Poisson brack¬ 
ets: 


(3.15) 


dW 

dt 


{H, W} . 
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We can also represent the dynamics via the proper time by using the repre¬ 
sentation dr = (I/ 7 )dt = (me 2 /H)dt, so that: 

dW dt dW H 


me“ 


{H, W} . 


dr dr dt 

Assuming a well-defined (invariant) rest energy (me 2 ) for the particle, we 
determine the canonical proper-time Hamiltonian K such that: 


{'<, = »'} , K | p=0 = H| p , 0 = 


me 2 . 


mcr 


Using 


[K, W} = 


' f a#' 

dW 

' H dH~ 

dw 

me 2 5p 

dx 

me 2 dx 

dp 


d_ 

dp 


H 2 

2 me 2 


+ ci 


dW d 


dx dx 


H 2 

2 me 2 


+ a 


dW 

dF’ 


we get that a = a' = \me 2 , so that (assuming no explicit time dependence) 


H 2 me 2 


K = and = {K, W} . 

2 dr 


2 me 2 


Since t is invariant during interaction (minimal coupling), we make the 
natural assumption that (the form of) K also remains invariant. Thus, if 
V c 2 p 2 + m 2 c 4 —>• y/c 2 Tv 2 + m 2 c 4 + V, where A a vector potential, V is a 
potential energy term and tt = p — -A. In this case, K becomes: 


7T 


K ~2^ + mC * + 2 me 2 


V 2 V sfc 2 ! r^”+rn?c^ 


+ 


me“ 


If we set Hq = \/c 2 7 t 2 + m 2 c 4 , use standard vector identities with Hq = 
meb — V, V x 7r = — |B, and compute Hamilton’s equations, we get: 


dx 

\ 1+ n 

7T 

mbc 

dr 

. #o_ 

m 

mbc — V 


7T U 

— => 7T = mu-u 

m 6 c 
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and 


(3.16) 


dp 

dr 


(■7T • V) 7T + |7T X B 


m 


+ 

t-H 

i _i 

- VU 

H 0 

m.c 2 

f 1 + ^l 

L Hoi 

vv- 

c 

mcb 




Further reduction, using the definition of E, with V = e<&, we have: 

V “ 


(3.17) 


dp e dA 
dr c dr 


G d A. 

x B — eV<F 

b or 0 


1 + 


mcb 


= eE + |u x B — eVd> 


V 

mcb 


It is clear that the additional term in equation (3.7) acts to oppose the force 
imposed by the charged particle part of the E held (i.e., — W). In order 
to see the physical meaning of the term, assume an interaction between a 
proton and an electron, where A = 0 and V is the Coulomb interaction, so 
that (3.7) becomes: 


(3.18) 


= -W- VV-^-r. 
b dr mcb 


Using Hq ~ me 2 , we see that lim r _>. ro u = 0 and lim r _>. ro a = 0, so that: 


(3.19) 


o = -VU- vu- 


V 


me* 


and the classical electron radius, ro, is a critical point (i.e., — VU — 
VF(h/mc 2 ) = 0). Thus, for 0 < r < ro, the force becomes repulsive. 
We interpret this as a fixed region of repulsion, so that the singularity r = 0 
is impossible to reach at the classical level. The neglected terms are attrac¬ 
tive but of lower order. This makes the critical point less than ro- Thus, in 
general, the electron experiences a strongly repulsive force when it gets too 
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close to the proton. This means that the classical principle of impenetra¬ 
bility, namely that no two particles can occupy the same space at the same 
time occurs naturally. It is this additional term that leads us to suspect that 
the electron may not act like a point particle in the s-states of hydrogen, 
where it has a finite probability of being at the center of the proton. 

The above observation also implies that, two electrons will experience an 
attraction if they can come close enough together as for example, at very 
low energies (temperatures). 

The Lagrangian representation reveals the close relationship to the non- 
relativistic case. If we solve for p, we get 


Vu e 

p = mu--—|—A. 

cb c 


Using this in K along with b 2 = u 2 + c 2 , we have 


g-t"" .-^) 2 + m^ + A + v( " M *.- y ) 


2 m 


2 me 2 


me* 


, 2 Vu 2 V 2 u 2 2 U 2 

= -mu~ -;-b -—rrv + rnc — 


be 2 mb 2 c 2 

From Cdr = p • fix — Kdr. we can write L as 

Uu 


_ Vb 

2 me 2 c 


C = 


mu — 


u + - A • u 


be 

Vu 2 V 2 u 2 


- ( U “ 2 - — + mv- 


+ me + 


Vb 


U 2 | 

2me 2 J 


i 2 e a 2 Vb V 2 

= -^mu + - A • u — me -b 


2 me 2 


1 - 


u 


b 2 


From this representation, it is clear that the neglect of second order terms 


gives us the non-relativistic theory. 
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4. Relativistic Quantum Theory 


The Klein-Gordon and Dirac equations were first discovered in early at¬ 
tempts to make quantum mechanics compatible with the Minkowski for¬ 
mulation of special theory of relativity. Both were partially successful but 
could no longer be interpreted as particle equations and a complete theory 
required quantum fields and the associated problems. For a recent discussion 
of other problems, one can consult [9j (see also M). 

In this section we introduce the canonical extension of the Dirac and 
square-root equations. Let A(x, t) and R(x) be given vector and scalar 
potentials and, after adding V (x) and making the transformation p —>• n = 
p-fA. 

To quantize our theory, we follow the standard procedure leading to the 
equation: 


3<F 

ih— = K<& = 

OT 


H 2 me 2 
2 me 2 2 


However, in addition to the Dirac Hamiltonian, there are two other possible 
Hamiltonians, depending on the way the potential appears with the square- 
root operator: 

/3-\/c 2 7t 2 — echTi ■ B + m 2 c 4 + V 


and 


c 2 7T 2 — echTi ■ B + (me 2 + /3V) 2 
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We have identified three possible canonical proper-time particle equations 
for spin-1 particles. (We also note that, these equations can be modified to 
apply to particles of any spin, by a minor change in the /3 matrix.) 


(1) The canonical proper-time version of the Dirac equation: 


(4.1) 


.M (tv 2 2 ehS-B 

ih— = \ - -h f3V + me 2 --- 

or 2m, 2mc 


+ 


Vat ■ tv ihot ■ W V 2 \ 


me 


2me 2 m.c 2 j 


T. 


(2) The canonical proper-time version of the square-root equation, using 
the first possibility: 


d'h f tv 2 ehS-B 
in—— = 


(4.2) 


dr 1 2m 2 me 


+ me + 


V 2 } 
2m,c 2 J 




V/3\Zc^ri^~^c^^^~+m 2 c^ /3y/c^Tv^~^cES^^~T^mJc^ 

-|— -\]> 4 - - V'fy. 

2 me 2 2m,c 2 

(3) The canonical proper-time version of the square-root equation, using 

the second possibility: 


(4.3) 


9T tv 2 9 

ih— = - - \- /3V + me - 

or 2m 


ehS • B V 2 

+ 


2m,c 2 me 2 


If V = 0, all equations reduce to: 


dT f tv 2 2 

= \^ + mc 


ehYi ■ B 
2 me 


T. 


The close relationship to the Schrodinger operator, makes it easy to see that, 
in all cases, K is positive definite. In mathematical terms, the lower order 
terms are relatively bounded with respect to tv 2 /2m. It follows that, unlike 
the Dirac and Klein-Gordon approach, we can interpret (4.1)-(4.3) as true 


particle equations. In the above equations, we have assumed that V is time 
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independent. (However, since A(x, t) can have general time-dependence, 
Vc 2 7r 2 — echT, ■ B + m 2 c 4 need not be related to the Dirac operator by a 
Foldy-Wouthuysen type transformation.) 

We plan to investigate the last two equations at a later time. In the next 
section, we focus on the canonical proper-time Dirac extension. 

4.1. The Dirac Theory. Since the Dirac equation forms the basis for QED, 
an important test of our proper-time extension is how well it compares to the 
Dirac equation in its description of the hydrogen spectrum. In this section, 
we compare the Dirac equation with the canonical proper-time extension for 
the Hydrogen atom problem. 

2 

If we let A = 0, Vo = —^r and consider the standard Dirac Hydrogen 
atom eigenvalue problem, 

A n ^n = H D 'Jf n , 


where A n is the n-th eigenvalue and 'V r) is the corresponding eigenfunction. 
For this case, if j is the total angular momentum and a is the fine structure 
constant, we have 

- 1/2 


A„ = mc~ 


1 + 


a 


n 


~\j + \\ + y(j + ¥) ~ a 2 


For the proper-time extension, with the same eigenfunction, we have 


E 4 / = 

^ n 


H 


D 


2 me 2 


+ 


mcr 


'F„ = 


A, 2 


2 me 2 


+ 


me 


'E, 
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In order to compare the two expressions to order 0(a 8 ), we begin with the 
following approximations: 

(1) (1 + x) _1 ~ 1 + x + x 2 + 0(x 3 ) 

(2) (1 + x)~ 1/2 ~ 1 — \x + |x 2 + 0{x 3 ) 

(3) (1 — x) 1/2 ~ 1 — ^x — |x 2 + 0(x 3 ). 


Let k = j + then our proper time extension can be written as: 



Using (2), we have: 
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and 
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With the last result, we now have: 
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It is now easy to see that, to order a 4 , \ n — E n = so that the E n values 
are systematically lower than the A,,, values. 

Table 1 below provides a relative comparison between the Dirac and 
proper-time extension compared with the experimental data for s-states, 
compiled by National Institute for Standards and Technology (NIST) of the 
US government. 


Table 1: Comparison with NIST data for s-states 


State 

Dirac 

Proper-time 

Nist 

A-DNIST 

A-PTNIST 

2s 

10.20439429 

10.20422448 

10.19881008 

.00558421 

.00541440 

3s 

12.09411035 

12.09393146 

12.08749443 

.00661592 

.00643603 

4s 

12.75550914 

12.75532871 

12.74853244 

.00697670 

.00679627 

5s 

13.06164150 

13.06146066 

13.05449789 

.00714361 

.00696277 
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As can be seen from the last two columns, the proper-time extension 
consistently provides results that a closer to the experimental data for all 
cases. (We could not compare the two ls-states with experiment, because 
of the NIST normalization for this state.) 


In Table 2, we see the same comparative results for the p, d and f-states. 

Table 2: Comparison with NIST data for p, d and f-states 


State 

Dirac 

Proper-time 

Nist 

A-DNIST 

A-PTNIST 

2p 0=1/2) 

10.20439429 

10.20422448 

10.19880553 

0.005588760 

0.005418952 

2p 0=3/2) 

10.20443957 

10.20426976 

10.19885089 

0.005588681 

0.005418870 

3p 0=1/2) 

12.09411035 

12.09393146 

12.08749292 

0.006617431 

0.006438537 

3p 0=3/2) 

12.09412377 

12.09394488 

12.08750636 

0.006617407 

0.006438512 

3d 0=3/2) 

12.09412377 

12.09394488 

12.08750634 

0.006617430 

0.006438535 

3d 0=5/2) 

12.09412824 

12.09394935 

12.08751082 

0.006617422 

0.006438528 

4p 0=1/2) 

12.75550914 

12.75532871 

12.74853167 

0.006977467 

0.006797044 

4f 0=7/2) 

12.75551763 

12.75533720 

12.74854038 

0.006976250 

0.006796820 


Thus, in all cases, the canonical proper-time extension of the Dirac equa¬ 
tion provides a closer approximation to the known experimental data for the 
Hydrogen spectra compared to the Dirac equation. In all cases, the changes 
are in the forth decimal place. This is insufficient to account for either the 
Lamb shift or the anomalous magnetic moment. 

2 

4.2. Future Direction. In what follows, let V = Vo = —p~. Based on 
our analysis of the square-root operator in the first section and the Dirac 
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operator in the second section, we are in the process of investigating the 
possibility that in s-states, the potential energy takes on the form: 


V = - 


(4.7) 
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where ro = e 2 /Me 2 . There are three possible choices for M: 

(1) The electron cannot be treated as a point particle in s-states of 
hydrogen, so that M = m, the mass of the electron and ro is the 
classical electron radius. 

(2) Neither the electron nor the proton can be treated as point particles 
in s-states of hydrogen, so that M = //, the reduced mass and ro is 
the classical mixed reduced radius. 

(3) The electron can be treated as a point particle in s-states of hydro¬ 
gen, but the proton cannot so that M = m p , the mass of the proton 
and ro is the classical reduced proton radius. 

It is clear that, at the zero-th order, we recover the Coulomb potential and 
the second term in (4.7) is a first order approximation. Assuming the first 
case, our eigenvalue problem becomes: 


(4.8) 
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As a first try, we set V = Vo for the terms containing a , use our first order 
approximation in the second term and V itself in the last term, so that 


r Zmc r 2m 


For the last term, we use the approximation: 

V 2 = 1 f A eKL- B | tt 2 Vo 
2me 2 2 | me 2 V m 2 c 3 m 2 c 2 me 2 V 
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2 1 [mc 2 J [ m 2 c 3 m 2 c 2 J J me 2 m 2 c 2 
Using these terms, we have: 

= + fiU>+ [i + «!5_il A EB 

\ P r r 2 j 2m [ P r r 2 \ 2 

nTr Vo a • 7r T iha- VVo T U n 2 _ ro2nft 
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= E 0 V + K'V, 
where 


K '= r o R r °] ^ [ r o r 0 l eltS-B r 0 M 2 

^.2 ” r 2m. r 2 P r 2me r m 2 c 2 


In conclusion, if this approach is as successful as we believe, we will still 
need to justify our approximation methods. 
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